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We study the eigenvalues of Dirac operators in QCD with two mass degenerate dynamical 
fermions. The gauge configurations have been obtained with HMC and the so-called Chirally 
Improved fermionic action. We compare eigenvalues obtained for the overlap Dirac operator on 
these configurations with those for the Chirally Improved (CI) operator (studied earlier). Results 
of Random Matrix Theory allow us to determine the chiral condensate. 
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1. Dirac operators: £> W i, D y, and D C i 

In this work we compare the properties of the eigenvalue spectra of three lattice Dirac oper- 
ators, namely the Wilson operator Dwi, the overlap operator Dov and the Chirally Improved (CI) 
operator, denoted as Dqi- The well known Dwi (for vanishing mass parameter m) is given by 

4 1 ±4 

D m (m,n) = £ (l-y^Z/^n) , (1.1) 

a la \i=±\ 

using the notation = — y^. The operator violates chiral symmetry, and thus does not satisfy the 
Ginsparg-Wilson equation. The second fermionic action studied here is defined through Dov [01 : 

D v = ^(l-y 5 sign(y 5 A)) , (1.2) 

with A = 1 s — a Dwi and < s < 2. For our calculations, we chose s = 1.8. 
This operator is an exact solution of the Ginsparg-Wilson equation 

y 5 D + Dy 5 =Dy 5 D , (1.3) 

and therefore implements chiral symmetry on the lattice. As a consequence of this the spectrum of 
Dqv lies exactly on a circle, the so-called Ginsparg-Wilson circle. 



The third Dirac operator we used, Dei, represents an approximate solution to (1.3). It is 
defined [||] as a truncated expansion of a most general solution of the Ginsparg-Wilson equation 
into 'paths' on the lattice of varying length. Taking paths up to infinite length results in an exact 
solution. Using this technique we can combine lower computer cost with - approximate, but good 
- chiral properties [ph. 



2. Dynamical Chirally Improved fermions 

For the following analysis of the spectral properties of the Dirac operators we use gauge fields 
with 2 dynamical flavors of fermions with degenerate masses. These were constructed with an 
HMC-algorithm, implemented with the Liischer-Weisz gauge action and the already mentioned CI 
fermionic action. More details can be found in [|t]. 

Due to the 'almost chiral' properties of this action, our HMC produces gauge-configurations 
which frequently tunnel between different topological sectors within one Markov chain. Table |] 
gives a short summary of the parameters of the gauge fields used in our analysis. 



3. Comparison of the spectral properties 

Inspection of the spectrum of a Dirac operator is a very good method to see how well chiral 
symmetry is approximated. In Fig. [j], f° r example, it can be clearly seen that the CI spectrum 
(right) deviates less from the ideal Ginsparg-Wilson circle than that of Dwi (left). 

Another information the spectrum provides is the topological charge Q top of a gauge configu- 
ration. According to the Atiyah-Singer index theorem it is possible to determine Q top by counting 
the zero modes of a Ginsparg-Wilson type Dirac operator according to their chirality, 

Q t0 p=n--n + (3.1) 
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run 


Z 3 x T 


Pi 


am 


a[fm] 


aniAw/ 


#confs. 


a 


12 3 x 24 


5.2 


0.02 


0.115(6) 


0.025 


73 


b 


12 3 x 24 


5.2 


0.03 


0.125(6) 


0.037 


52 


c 


12 3 x 24 


5.3 


0.04 


0.120(4) 


0.037 


55 


d 


12 3 x 24 


5.3 


0.05 


0.129(1) 


0.050 


40 



Table 1: The parameters for the runs used for this work. L 3 x T denotes the extent of the lattice in units of 
the lattice spacing a, am the bare mass parameter of Dq, am^wi the quark mass calculated via the axial Ward 
identities, and j3i the first gauge coupling of the set of three LW-couplings used The configurations are 
separated by 10 HMS-trajectories. 
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Figure 1: Left: Part of the spectrum of Dwr for configuration nr. 100 of run a. Right: The 50 smallest 
eigenvalues of Dei for the same gauge-configuration. 



with n± denoting the number of eigenmodes |v) with chirality (v\y$ \ v) = ±1 corresponding to 
eigenvalues X = 0. When treating Dqi, due to its approximate nature, these eigenvalues are not 
exactly zero, but scatter on the real axis. For the CI operator (and the Wilson operator) one may 
relate the number of real modes to the topological sector. We determine Q top by setting n± as 
the number of eigenmodes \w) corresponding to real eigenvalues with chirality (w|y5|w) ^ 0, 
respectively. (For eigenvalues not on the real axis we numerically find (w | 75 \w ) = as expected.) 
When comparing the topological charge of a configuration calculated with the exactly chiral Dov 
and the approximately chiral Dei, and we find approximate agreement (cf. Fig. fy. The differences 
mostly originate from missed eigenvalues far inside the Ginsparg-Wilson circle, which are not 
recovered in our method of calculating eigenvalues, namely the program package ARPACK; this 
tool only computes the eigenvalues with the smallest absolute value with respect to a defined origin. 
When computing eigenmodes of the overlap operator, the situation is different. There the sector 
depends on the value of s added to the diagonal part of the kernel operator (Wilson in our case), 
which has frequent inner real modes that are missed in that overlap projection, depending on s. 

If one does not take into account any normalization and directly compares eigenvalues of Dq\ 
with those of Dov as defined in (1.2), an interesting behavior can be seen. In cases where not only 
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Figure 2: Qrop.oy vs. Q t0 p.ci for the gauge configurations of run a. The numbers state how many configura- 
tions show a particular combination of topological charges. 



the topological charge Q top is identical for both operators but also the number of real eigenvalues 
no, the first few eigenvalues are in direct correspondence (cf. Fig. |3[ left). However, if no > Q top 
for the CI operator, i.e., if eigenmodes with opposite chirality cancel each other with respect to 
Qtop, this correspondence is lost (cf. Fig. ||, right). When put into the overlap operator, the real 
modes seem to "move" up resp. down the imaginary axis (although not all of them close to the real 
axis), thereby increasing the eigenvalue density near the origin. For a similar observation cf. [Q]. 
This enhancement leads to a higher value for the (bare) chiral condensate when calculated directly 
with this definition of Dqv- 
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Figure 3: Eigenvalues Xpv of £>ov an d eigenvalues X^i of Dq m me complex plane. Left: configuration nr. 
1 10 (run (a)): Q top — —2, n + —2 for both operators Right: configuration nr. 100 (run (a)): Q top = —2, but 
n+,ci = 3, n-a = 1. 
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Figure 4: The average value ipr(X) over the complex plane for the gauge configurations of run (a). 

4. Localization of eigenmodes 

For the CI operator, we also calculated the inverse participation ratios (ipr) of the eigenmodes 
for small eigenvalues, given by 

ipr(X)=V Wj\(*, a,c)*v(*,a,c) I . (4.1) 

x \a,c I 

The inner sum is over the color indices c as well as the Dirac indices a, the outer sum over the 
space-time indices x. V denotes the space-time volume in lattice units, and v(x, a,c) the eigenmode 
of Dei corresponding to the eigenvalue X. This quantity is a good measure for the localization 
properties of one eigenmode, with ipr(X) = 1 for a non-localized mode and ipr(X) = V for a mode 
concentrated on only one lattice point. To compute a suitable average over one HMC-run, we 
divided the complex plane into a grid and calculated the average ipr as 

Jp7=- £ i pr (X), (4.2) 

with A^ being one square in the complex plane and n the number of eigenvalues in A^ . 

In Fig. ^| we see, as expected, that ipr increases along the real axis, and in general is higher 
for eigenvalues inside the Ginsparg-Wilson circle (cf. also [Q]). For Dov one expects ipr to be 
symmetric with respect to the real axis, but this is not the case for Dq\ since it is not a normal 
matrix operator. 

5. Comparison with random matrix theory 

As the overlap operator implements chiral symmetry on the lattice exactly, the distribution of 
the smallest eigenvalues (in leading order ChPT) is analytically given by the well known results of 
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random matrix theory (RMT) for the chiral gaussian unitary ensemble [Q], at least in the £-regime 
in the microscopic limit. Following the procedure in [Q], where these calculations were done on 
the same configurations for Dq\, we thus compare our distributions for Dov to RMT. 

Jp(AEV) Jp(ASV) 




12345 2468 10 



Figure 5: The cumulative distribution given by RMT compared to the distribution of the smallest eigenvalues 
of the overlap operator. Left: Q top = 0, and k=l (smallest eigenvalue); Right: Q top = 0, and k=2 (second 
smallest eigenvalue), both for run (a). 

A similar analysis for dynamical overlap configurations was done in [|| ^ and for 2-flavor 
staggered configurations in JTOj]. For the fits that determine the chiral condensate E, we use the 
cumulative eigenvalue distributions and look at the smallest and second smallest eigenvalue in the 
topological sectors Q top = and \Q t0 p \ = 1- 

The values for £ for run (a)-(d), determined from the spectra of Dov, are given in Table ||. 
The fit was done using the Kolmogorov-Smirnov test, the errors computed via statistical bootstrap. 

6. Problems and issues 

The resulting values for the condensate still have to be renormalized. This can be done by 
determining the renormalization constant Z$ by standard tools of non-perturbative renormalization. 
In the case of Dq\ the renormalization constants for the dynamical case have been computed in 
recent work JTl|]. For Dov, the situation is more involved, and the renormalization depends on the 
value of s. The weak coupling expansion for small momenta p has a behavior Dov = i YuPn/s + 
0(p 2 ). This changes for the interacting case in a non-trivial way. Comparing the results for the 
bare condensate for CI [§] with overlap, we expect Zs t ov = Zs,ci^ci/^ov ~ 0.64Zsc/- 

The dependence on the physical quark mass is also worth being explored. To this point we 
assume that the AWI-mass computed for our dynamical gauge fields is the same as the quark mass 
entering the RMT-fits. In [ |l0| ] the overlap quark mass has been determined using the distribution 
of topological charge, and found to be different from the value calculated for the underlying gauge 
fields. We do not have sufficiently high statistics to follow this approach, but nevertheless the sea 
quark mass of the overlap operator should be adjusted, such that the pion mass (or the AWI-mass) 
computed with the overlap operator agree with the pion mass of the CI operator. 
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HMCrun 


k 




#confs. 


-(E) 1 / 3 MeV 1 / 3 


a 


1 





12 


338(6) 




1 


1 


25 


332(3) 




2 





12 


310(3) 




2 


1 


25 


319(3) 


b 


1 





7 


353(10) 




1 


1 


8 


350(9) 




2 





7 


362(5) 




2 


1 


8 


330(6) 


c 


1 





17 


350(6) 




1 


1 


12 


346(9) 




2 





17 


340(3) 




2 


1 


12 


322(5) 


d 


1 





5 


365(18) 




1 


1 


11 


370(10) 




2 





5 


348(14) 




2 


1 


11 


346(2) 



Table 2: Results for the value of the bare condensate E, for all runs of the HMC. The value of k refers to the 
smallest (1) or second smallest (2) (imaginary part of the) eigenvalue. 
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